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Abstract. The main objective of this work is to extend the averaging method 
for studying the periodic orbits of a class of differential equations with discon- 
tinuous second member. Thus, overall results are presented to ensure the 
existence of limit cycles of such systems. Certainly these results represent new 
insights in averaging, in particular its relation with non smooth dynamical 
systems theory. An application is presented in careful detail. 

1. Introduction and statement of the main results 

The discontinuous differential systems, i.e. differential equations with discon- 
tinuous right-hand sides, is a subject that has been developed very fast these last 
years. It has become certainly one of the common frontiers between Mathematics, 
Physics and Engineering. Thus certain phenomena in control systems, impact in 
mechanical systems and nonlinear oscillations are the main sources of motivation 
of their study, see for more details Teixeira [21] . 

The knowledge of the existence or not of periodic solutions is very important 
for understanding the dynamics of the differential systems. One of good tools for 
study the periodic solutions is the averaging theory, see for instance the books of 
Sanders and Verhulst [19] and Verhulst [23]. We point out that the method of 
averaging is a classical and matured tool that provides a useful means to study the 
behaviour of nonlinear smooth dynamical systems. The method of averaging has 
a long history that starts with the classical works of Lagrange and Laplace who 
provided an intuitive justification of the process. The first formalization of this 
procedure was given by Fatou in 1928 [8]. Very important practical and theoretical 
contributions in the averaging theory were made by Krylov and Bogoliubov [2] in 
the 1930s and Bogoliubov [T] in 1945. The principle of averaging has been extended 
in many directions for both finite- and infinite-dimensional differcntiablc systems. 
The classical results for studying the periodic orbits of differential systems need at 
least that those systems be of class C 2 . Recently Buica and Llibre [4] extended 
the averaging theory for studying periodic orbits to continuous differential systems 
using mainly the Brouwer degree. 

The main objective of this paper is to extend the averaging theory for studying 
periodic orbits to discontinuous differential systems using again the Brouwer degree. 

Let D be an open subset of K™. We shall denote the points of R x D as (t,x), 
and we shall call the variable t as the time. Let /i:Rx£>— ^RbeaC 1 function 
having the £ 1 as a regular value, and let M = /i _1 (0). 
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Let X, Y : R x D — > W 1 be two continuous vector fields. Assume that the func- 
tions h, X and Y are T-periodic in the variable t. Now we define a discontinuous 
piecewise differential system 



(1) x'{t) = Z(t,x) 



X(t,x) if h(t,x)>0, 
Y{t,x) if h(t,x)<0. 



Denote Z = (X,Y). 

Here we deal with a different formulation for the discontinuous differential system 
([l]). Let sign(u) be the sign function defined in R \ {0} as 



sign(w) 



1 if u > 0, 
-1 if u<0. 



Then the discontinuous differential system ([T]) can be written using the function 
sign(u) as 

(2) x'(t) = Z(t,x) = F 1 (t,x)+sign(h(t,x))F 2 (t,x), 

where 

F 1 (t,x) = ^(X(t,x)+Y(t t x)) and F 2 (t, x) = - (X(t, x) - Y(t, x)) . 

To work with the discontinuous differential system ^ we should introduce 
the regularization process, where the discontinuous vector field Z(t,x) is approx- 
imated by an one-parameter family of continuous vector fields Zs(t,x) such that 
lim^o Z§ = Z(t,x). 

In |20j Sotomayor and Teixeira introduced a regularization for the discontinuous 
vector fields in R 2 having a line of discontinuity and, using this technique, they 
proved generically that its regularization provides the same extension of the orbits 
through the line of discontinuity that the one given by the Filippov's rules, see [8]. 
Later on Llibre and Teixeira |18j studied the regularization of generic discontinu- 
ous vector fields in R 3 having a surface of discontinuity, and proved that lim^^o %S 
essentially agrees with Filippov's convention in dimension three. Finally, in |21j 
Teixeira generalized the regularization procedure to finite dimensional discontinu- 
ous vector fields. 

In |16j Llibre, da Silva and Teixeira studied singular perturbations problems in 
dimension three which are approximations of discontinuous vector fields proving 
that the regularization process developed in [TS] produces a singular problem for 
which the discontinuous set is a center manifold, moreover, they proved that the 
definition of sliding vector field coincides with the reduced problem of the corre- 
sponding singular problem for a class of vector fields. 

In all these regularizations (except the one using singular perturbation theory) 
a transition function is used to average the vector fields X and Y on the set of 
discontinuity in order to get a family of continuous vector fields that approximates 
the discontinuous one. 

A continuous function <j) : R — > M is a transition function if (j){u) = —1 for 
u < —1, 4>(u) = 1 for x > 1 and <fr'(u) > if u e (—1, 1). The (/>-regularization of 
Z = (X, Y) is the one-parameter family of continuous functions Zg with S £ (0, 1] 
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given by 

Z s (t, x) = \ (X(t, x) + Y(t, x)) + l -h(Ht, x)) (X(t, x) - Y(t, x)) , 

with 

(3) &(«) = *(|). 

Note that for all (t, x) G (K x D)\M we have that lim^o Zs(t, x) = Z(t, x). 

The formulation © of the discontinuous differential system ([1]) admits a natural 
rcgularization. Define the transition function <f> as 

{1 if «>1, 
u if - 1 < u < 1, 
-1 if u<-l. 

Let cj>s : K — > M be the continuous function defined in (j3J). It is clear that 

(5) lim 4>s(u) = sign(u), 

and 

Z s {t,z) = Fi(t,x) + <j> s {h{t,x))F 2 (t,x), 
is the 4>-regularization of the discontinuous differential system ([1}. 

If S = {0} x D, then doing a time translation (if necessary) we always can 
assume that £ ^ M.. As usual V/i denotes the gradient of the function h, and V x /i 
denotes the gradient of the function h restricted to the variable x. 

Our main result on the periodic orbits of the discontinuous differential system [2] 
using averaging is given in the next theorem. Its proof uses the theory of Brouwer 
degree for finite dimensional spaces (see the appendix A for a definition of the 
Brouwer degree ds(/, V, 0)), and is based on the averaging theory for non-smooth 
differential system stated by Buica and Llibre [I] (see Appendix B). 

Theorem 1. We consider the following discontinuous differential system 

(6) x'(t) =eF(t,x) +e 2 R{t,x,e), 
with 

F(t,x) = Fi{t,x) + sign(h(t,x))F 2 (t,x), 

R(t, x, e) — Ri(t, x, e) + sign(h(t, x))R 2 (t, x,e), 

where F 1: F 2 : M. X D -» W 1 , R X ,R 2 :lx D x (-e , e ) -> R" and h : K x D -> R 
are continuous functions, T -periodic in the variable t and D is an open subset of 
K™. We also suppose that h is a C 1 function having as a regular value. Denote 
by M = /i -1 (0), by £ = {0} x D <£ M, by E = ^ 0. and its elements by 

z = (0,2!) g X. 

Define the averaged function f : D — > M" as 

(7) /(x) = / F(t,a;)dt. 

Jo 

FFe assume the following three conditions. 

(i) Fi, F 2 , Ri, R 2 and h are locally L-Lipschitz with respect to x; 

(ii) for a g So with /(a) = ; there exist a neighbourhood V of a such that 
f(z) ± for all zeV\{a} and d B (f,V,0) + 0. 
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(iii) If dh/dt ^ 0, then for all (i, z) 6 Ai we have that (dh/dt)(t, z) ^ 0; and if 
dh/dt = 0, then (V x h,Fi) 2 - (V x h,F 2 ) 2 > for all (t,z) G [0, T] x M. 

Then, for \e\ > sufficiently small, there exists a T~periodic solution x(-,e) of 
system ([6]) such that x(0,e) — > a as e — > 0. 

Theorem Q] is proved in section [21 

We remark that condition (iii) of Theorem [T] changes according with the fact 
that the hypcrsurfacc of discontinuity M. depends or not on the time . 

The averaging theory has already been used in discontinuous dynamical systems. 
Thus, for instance, Llibre, Novaes and Teixeira [13] have used the averaging theory 
to provide sufficient conditions for the existence of periodic solutions of the planar 
double pendulum with discontinuous perturbation. Llibre and Teixeira [17j have 
used the averaging theory to provide lower bounds for the maximum number of limit 
cycles for the m-piecewise discontinuous polynomial differential equations. In [5] 
Cardin, Carvalho and Llibre have used the averaging theory to study the bifurcation 
of limit cycles from the periodic orbits of two and four dimensional linear center in 
R n perturbed inside in a class of discontinuous pieccwise linear differential systems. 
Finally, in [15] Llibre and Rong have used the averaging theory to study the number 
of limit cycles of the discontinuous piecewise linear differential systems in K 2n with 
two zones separated by a hyperplane. 

In Theorem [T] we have extended to general discontinuous differential systems the 
ideas used in the previous mentioned papers for particular discontinuous differential 
systems. 

Now we shall do an application of Theorem [T] to the discontinuous piecewise 
linear differential systems separated by a straight line. These differential systems 
have been studied recently by Han and Zhang [TO], and Huan and Yuang [TTj . 
among other papers. In [10] some results about the existence of two limit cycles 
appeared, so that the authors conjectured that the maximum number of limit cycles 
for this class of piecewise linear differential systems is exactly two. This conjecture 
is analogous to Conjecture 1 in the discussion of Tonnelicr in [22] . However, by 
considering a specific family of discontinuous PWL differential systems with two 
linear zones sharing the equilibrium position, in |11] strong numerical evidence 
about the existence of three limit cycles was obtained, and a proof was provided by 
Llibre and Ponce [14]. This example represents up to now the first discontinuous 
piecewise linear differential system with two zones and 3 limit cycles surrounding 
a unique equilibrium. Now we shall provide a new proof of the existence of these 
three limit cycles using the averaging theory through our Theorem [JJ 

In polar coordinates (r, 9) given by x = r cos 9 and y = r sin 9, the planar discon- 
tinuous piecewise linear differential system with two zones separated by a straight 
line corresponding to the system studied in the paper |14j is 



where we have multiplied the right hand side of the system of [14j by the small 
parameter e. Our main contributions in this application is to provide the explicit 



(8) 
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analytic equations defining the limit cycles of the discontinuous piecewise linear 
differential system with two zones (|8]). 

Theorem 2. Any limit cycle of the discontinuous piecewise linear differential sys- 
tem with two zones ((U) which intersects the straight line x = 1 in two points (rg, 9q) 
and (r±, 9\) with —tt/2 < Qq < < 9\ < ir/2, r^ cos 9k = 1 for k = 0, 1, and r<j > 1 
and 9\ must satisfy the following two equations 

/ 19(gi-go) \ a 1 n 
exp — r cos 9\ — 1 = 0, 



V 50 

1 + ^ arctan \—= sec 9 Q (23 cos 9 — 100 sin # 



50 5 \15 

(9) 1/1 

v ' —- arctan I — sec 6*i (23 cos 6*i — 100 sin 6*i 

V 

-ilog(|4623cos(26> ) + 2300 sin(26» ) - 5377|) 

1 27T 

+ -log(|4623cos(26ii) + 2300sin(26li) - 53771) = 0, 

2 5 

where 9q = arccos(l/ro) — tt, and the determination of the arctan is in the interval 
(— 7r/2,7r/2) and of the arccos in the interval (0,7r). 

On the other hand, any limit cycle of system ([8]) which intersects the straight 
line x = 1 in two points (ro,9o) and (r±,9i) with < 9q < 9\ < tt/2, r^cos^ = 1 
for k = 0, 1, and ro > 1 and 9\ must satisfy the equations but now in both 
equations 9q = arccos(l/ro). 

We recall that a limit cycle of system (j8]) is an isolated periodic orbit of that 
system in the set of all periodic orbits of the system. It is well known that the 
study of the limit cycles of the differential systems in dimension two is one of the 
main problems of the qualitative theory of differential systems in dimension two, 
see for instance the surveys of Ilyashenko [9] and Jibin Li [12] . 

In fact, as we shall see in the proof of Theorem [2] equations ([9]) have three 
solutions, providing the three limit cycles of Figure [TJ 



2. Proof of Theorem Q] 

For proving Theorem [1] we need some preliminary lemmas. In what follows we 
always assume that all hypotheses of Theorem [T] hold. 

Lemma 3. If dh/dt ^ 0, then for z e D and t e [0,T] the map h z : t i-> h(t,z) 
has finitely many zeros in [0, T] . 

Proof. From the definition of h we know that the map h z : t i-> h{t,z) is C 1 . 
Suppose that there exist an increasing sequence (ii)i S N C [0,T] such that h z (ti) = 
for all i 6 N. Then there exist a convergent subsequence (ti^j^ such that 
ti. — > i G [0,T]. By continuity of h z . it follows that h z (t) = 0, which implies that 
(i, z) € A4. We know, by the Mean Value Theorem, that for each j G N there exist 
Sj G (ti-,ti. +1 ) such that h' z (sj) = 0. Since Sj — > t, then h' z {t) = 0, which is a 
contradiction with hypothesis (iii). Hence, the map h z : t n- h{x, z), for t G [0,T], 
vanishes only in a finite subset T C [0, T]. □ 
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Figure 1 . The three limit cycles surrounding the origin. 

In the hypothesis (ii) of Theorem Q] we use the Brouwer degree of the averaged 
function /. We need to show that the function / is continuous in order that the 
Brouwer degree will be well defined, for more details see Appendix A. 

Lemma 4. The averaged function (J7J) is continuous in Y>q. 

Proof. Fix zq G So- Since the map h Zo : t h-> h(t, zq) is C 1 and (TSLxD)\M is an open 
subset ofRxD. If dh/dt ^ 0, by LemmaEl for t G [0, T], we know that h(t, z ) = 
only for a finite subset T = {h,t 2 , ■ ■ ■ ,t k -i} of [0,T]. If dh/dt = 0, thenT = 0. Set 
to = and tk = T. Now h(t, Zo)h(s, zq) > for t, s G (U, ii+i) and i = 0, 1, • • ■ , k—1, 
i.e. h(t, Zq) does not change its sign in (i,, U + i) for i = 0, 1, • • ■ , k — 1. Now, for 
z G D in some neighborhood of zq, we estimate 

|/C0-/(*>)| < / \F 1 (t,z )-F 1 (t t z)\dt 

|sign(/i(t, z ))F 2 (t, z ) - sign(/i(i, z))F 2 {t, z)\dt 



< TL\z -z\+[ F(t,z ,z)dt 

Jo 

< TL\z -z\+^2J F(t,z ,z)dt, 



where F(t,z ,z) = \siga(h(t, z ))F 2 (t, z ) - sign(/i(t, z))F 2 {t, z)\. Set = 

and /7 = [t; + 7, — 7] for i = 1, 2, • • • , k — 1. Note that there exists 70 > such 

that J/ C Ii for all 7 < 70- The set I? x {^o} C 5 is compact, then there exist a 
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ball B{zo,r 1 ) such that K] = x B(z , r 7 ) C S. We can assume that r 7 < 7. Set 
Ml = sup{#(t,z ,^) : (t,z) € i<7}. Note that 

F(t,z ,z) < |sign(/i(t,z ))F 2 (t,z )| + |sign(/ l (t,z))F 2 (t,z)| 
= \F 2 (t,z )\ + \F 2 (t,z)\, 



if we take M = sup{|F 2 (t, z )\ + \F 2 (t,z)\ : (t,z) e [0,T] x 5(z ,r 7 )}, then M? < 
M for alH = 1, 2, • • • , k - 1. Thus 

+, #(t,z ,z) = / F(t,z ,z)+ f ^ F(t,z ,z)+ f + F(t,z ,z) 

< TL\zo-z\ + 2yM? <(TL + 2M ft. 

Hence 

|/(2o) - /(z)| < TL7 + k(TL + 2M )7 < 0>7- 
Therefore, given e > we can choose 7 = e/Cq, then for all z £ B(zq, r T ) we have 
|/(2o ) — /( z )l < £• This implies that / is continuous in z € So- □ 

Let a be the point in hypothesis (ii) of Theorem [TJ By Lemma [4] there exists a 
neighborhood V of a such that / is continuous in V. Hence, by Theorems |9] and 
1101 (see Appendix A) , there exists a unique map that satisfies the properties of the 
Brouwer degree for the function f(z) with z E V, because ^ f(dV). This map is 
denoted by d B (f,V,0). 

Lemma 5. For \e\ sufficiently small the solutions of system ([6]), in the sense of 
Filippov, are uniquely defined in any compact set [0, T] x Di contained in [0, T] x D. 

To prove Lemma [5] we will need the following proposition, that has been proved 
in Corollary 1 of section 10 of chapter 1 of [8]. Define 

S + = {(t,x) eRxD: h(x,t) > 0}, 
S~ = {(t,x) elxfl: h(x,t) < 0}. 
Note that RxD = S~l)Mll S + . 

Proposition 6. For every point of the manifold A4 where (Xh)(Yh) > 0, there is 
a unique solution passing either from S~ into S + , or from S + into S~ . 

Proof of Lemma [21 System (O can be written as the autonomous system 

(X(r,x) if h(T,x)>0, 
~\Y(t,x) if h(r,x)<0, 



inRxD, where 
X(t,x) = 

Y(t,x) = 



1 

li £(F 1 (r, x) + F 2 (t, x)) + e 2 {R 1 {r, x, e) + R 2 (t, x,e)) j 

1 

s(Fi(r,x) - F 2 (t,x)) + e 2 (Ri(T,x,e) - R 2 (t,x,e)) 
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dh^ 2 



If dh/dt^ 0, then we have that 

(Xh)(Yh) = (Vh,X)(Vh,Y) = (t,x) + 0(e). 

By hypotheses (iii) we know that if (t,x) E M, then (dh/dt)(t,x) ^ 0. Hence, for 
|e| sufficiently small, we conclude that (Xh)(Yh)(t, x) > 0, for every (t,x) in any 
compact set [0,T] x D\ with D\ C D. Applying Proposition [5] the result follows. 

If we assume that dh/dt = 0. Then 

(Xh)(Yh) = (V/i, X)(Vh, Y) = e 2 ({W x h, F±) 2 - (V x h, F 2 ) 2 ) + 0(e 3 ). 

By hypotheses (iii) for |e| sufficiently small we conclude that (Xh)(Yh)(t, x) > 
for every (t,x) in any compact set [0,T] x D\ with D\ C D. Again applying 
Proposition [S] the lemma is proved. □ 

Instead of working with the discontinuous differential system ([6]) we shall work 
with the continuous differential system 

(10) x'(t) = eF s {t,x) + e 2 R 5 {t 1 x,e), 

where 

F s (t, x) = F\ (t, x) + 4> s {h{t, x))F 2 (t, x), 
R$(t,x,s) = Ri(t,x,e) + 4>s(h(t,x))R 2 (t,x,s), 

and 4>s : R — > R is the continuous function defined in ^ and and satisfying 
©• 

For system (|T0|) the averaged function is defined as 

fs(z) = / F s {t,x)dt. 



We need to guarantee that hypothesis (i) of Theorem [TT] (see appendix A) holds 
for the functions Fs and R$. For this purpose we prove the following lemma. 

Lemma 7. For 6 E (0,1] the function 4>s : R — > R defined in ^ with <f> given by 
((4]) is globally 1 / '<$ -Lipschitz; i.e. for allui,U2 E R we /iawe £/ia£ ^(iti) — </>,5(ii 2 )| < 
(l/*)|«i-u 2 |. 

Proo/. If Ul < -5 < 8< u 2 , then |<fo(ui) - <fo(u 2 )| = 2 = (l/<*)2(5 < (l/<J)|iti -u 2 |- 

If Ui,U2 < or ui, u 2 > (5, then |0a(iti) - <j>s{u 2 )\ — < (l/5)]wi - u 2 |. 

Assume that u\ E (— <5, (5). If \u 2 \ < S, then \(ps(ui) — 4>s(u 2 )\ = (l/5)\ui — u 2 \; 
and if |u 2 | > S, then <As( u 2)| < max{|l/5|, |l/u 2 |}|ni— u 2 \ < (1/S)\ui— u 2 \. 

This completes the proof of the lemma. □ 

Proposition 8. For 5 £ (0, 1] £/ie functions Fs and Rs are locally Lipschitz with 
respect to the variable x. 

Proof. Let K C D be a compact subset. Denote M = sup{|F 2 (t, x)| : (t, x) E 
[0, T] x A'}, which is well defined by continuity of the function (t,x) i— > |.F 2 (f,a;)| 
and compactness of the set [0, T] x K. For and x 2 in K where F\ and h are 
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locally L-Lipschitz and by Lemma [7j we have 



\F s (t,x 1 )-F s (t,x 2 )\ 



\F 1 (t,x 1 )-F 1 {t,x 2 ) 

+fo{h{t,x l ))F 2 {t 1 x 1 ) - <t> s {h{t,x 2 ))F 2 {t,x 2 )\ 
\F 1 {t,x 1 )-F 1 {t,x 2 )\ 

+\<j) 5 {h{t,x 1 ))F 2 {t,x l ) - <t) 5 {h{t,x 2 ))F 2 {t,x 2 )\ 
L\ Xl - Xl \ + |^ 4 (/i(t,a;i))||F2(t,ii) - F 2 {t,x 2 )\ 

+\F 2 (t, X2 )\\MHt,xi)) ~ MKt,x 2 ))\ 

M 

2L\x! - x 2 \ + —\h(t,xi) - h{t,x 2 )\ 



< 



< 



< 



< 



2L H — |xi - x 2 \ = L s \x\ - x 2 \ 




The proof for R$ is analogous. 



□ 



Now we are ready to prove Theorem [T] 

Proof of Theorem QJ We will study the Poincare maps for the discontinuous differ- 
ential system (JB|) and for the continuous differential system (JTUJ). For each z G D, 
let xs{t, z, e) denote the solution of system (flu]) such that 2:5(0, z, e) = z. 

Since the differential system (|10[) is T-periodic in the variable i, we can consider 
system (flO)) as a differential system defined on the generalized cylinder S 1 x D 
obtained by identifying E = {(t, x) : r = 0} with {(r, x) : r = T}, see Figure [21 On 
this cylinder E is a section for the flow. Moreover, if z E D is the coordinate of a 
point on E, then we consider the Poincare map P$(z) = xs(T, z, e) for the points z 
such that xs(T, z,e) is defined. 

Observe that there exists Eq > such that, whenever e € [— So,£o], the solution 
x$(t, z, e) is uniquely defined on the interval [0, T]. Indeed, if (t~,t+) is the maximal 
open interval for which the solution passing through (0,2) is defined. Now we 
shall apply the local existence and uniqueness theorem for the solutions of these 
differential, see for example Theorem 1.2.2 of [19]. Note that we can apply that 
theorem due to the result of Proposition [H Hence, by the local existence and 
uniqueness theorem we have that £+ > h z and h z = inf {T, d\m(e)} where m(e) > 
\eFs(t, x) + e 2 Rs(t, x, e)\ for all t e [0, T], for each x with \x — z\ < d and for every 
z G D. When |e| is sufficiently small, m(e) can be arbitrarily large, in such a way 
that h z = T for all z e D. Hence, for e G [— £o,£o]: the Poincare map of system 
(|T0|) is well defined and continuous for every z G D. 

From the definition of the Poincare map Pg(z) its fixed points correspond to 
periodic orbits of period T of the differential system (flU)) defined on the cylinder. 

We can define in a similar way the Poincare map P £ {z) of the discontinuous 
differential system ([6]) . More precisely, by Lemma [5] and by the continuous depen- 
dence of solution of the differential equation with discontinuous righthand side on 
the initial data (see for instance Corollary 1 of section 8 of chapter 2 of [8] ) , such 
Poincare map is continuous. Again the fixed points of P e {z) correspond to periodic 
orbits of the discontinuous differential system ([6j) . 

By definition the continuous differential system (|10|) is C 2 in the variable e. So 
we do the Taylor expansion of the Poincare map of system (|10|) around e up to 
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nT=0 
S 




FIGURE 2. Generalized cylinder. 

order two, and we get 

(11) P!(z) = z + ef s (z) + 0(e 2 ), 

where fs(z) is the averaged function of the continuous differential system ffTUf) . for 
more details see for instance [3]. 

The referred Poincare map is the composition of the Poincare maps of the con- 
tinuous differential systems, so it is continuous. Therefore, we can repeat the ar- 
guments of [1] and we again have that 

(12) pc(z) = z + e.f(z) + 0(e 2 ), 

where f(z) is the averaged function of the discontinuous differential system 

Due to {5} the pointwise limit of the Poincare map P$(z) of system (|10[) . when 
S —> is the Poincare map P £ (z) of system ([6]). Consequently, from (|TT|) and (|T2|) 
we obtain that the pointwise limit in E of the function f$, when S — > is the 
function /. 

Let a G V be the point satisfying hypotheses (ii) of Theorem [T] Therefore 
f(z) ^ for all z 6 V\{a}. Define /o = f\y, we know that fo is continuous by 
Lemma |H Then, we consider the continuous homotopy {fs\y,0 < <5 < 1}. We 
claim that there exists a <5o G (0, 1] such that ^ fs(dV) for all 6 G [0, So]. Now we 
shall prove the claim. 

As usual N denotes the set of positive integers. Suppose that there exists a 
sequence (z m )meN in dV such that fj_(z m ) = 0. As the sequence (z m ) is contained 
in the compact set dV, so there exists a subsequence (zm^gN such z me — > zq g dV. 
Consequently we obtain that f(zo) = 0, in contradiction with the hypotheses (ii) 
of Theorem [TJ Hence, the claim is proved. 

From the above claim and the property (iii) of Theorem [3] (see Appendix A) we 
conclude that rfs(/<5,V,0) ^ for < S < Sq. Therefore, by the property (i) of 
Theorem [H] we obtain that G fs(V), so there exists ag G V such that fs(as) = 0. 
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Since, by continuity, there exists the lim^o as and it is a zero of the function 
fo = f\v ■ This limit is the point a of the hypotheses (ii) of Theorem [TJ because a 
is the unique zero of fo in V. 

In summary, in order that for every S G (0, 5q] the averaged function fs satisfy 
the assumptions (ii) of Theorem ITT1 (sec Appendix B). So it only remains to show 
that in V we have that fs(z) ^ for all z G V \ {as}. But this can be achieved in 
a complete similar way as we proved the above claim. Hence, by Proposition [8] for 
every 6 G (0, So) the continuous differential system (flu]) satisfies all the assumptions 
of Theorem ITTl Hence, for |e| sufficiently small there exists a periodic solution 
xg(t,e) of the continuous differential system ()10|) such that zrs,e) := xg(0,e) — > as 
when e — > 0. 

Now, from pip the point z^ :E ) is a fixed point of the Poincare map Pg(z), i.e. 
-P|( z (<5,£)) = z (S,e)- Since lim.g-,.0 Pg (z) = P £ (z), it follows that z e = lim^o Zr$ je ) is 
a fixed point of the Poincare map P e (z). So, the discontinuous differential system 
© has a periodic solution x(t,e) such that z e = x(0,e) — > a as e —> 0. Therefore 
the theorem is proved. □ 



In this section we shall prove Theorem [2 by applying Theorem [T] to the discon- 
tinuous differential system (j8|). So, we must compute the integral (0, which for 
system ((8]) becomes 



where the function F(6, r) is given in (|8|) . 

The solution of the differential system ([8]) in the half-plane x = rcos# > 1 
starting at the point (ro,#o) with rocos#o = 1 and 9q £ (— 7r/2,0) is 



Therefore, at the point (ri,#i) with ricos^i = 1 and B\ G (0,7r/2) we have that 



This equation coincides with the first equation of ([5]). 

Now computing the integral (fT5)) we obtain exactly the right hand side of the 
second equation of © multiplied by r. According to Theorem [T] we must find the 
zeros of this last expression. Since r cannot be zero the equation for the zeros is 
reduced exactly to the second equation of ©. In short, by Theorem [1] we have 
proved that a periodic orbit of system ((5J intersects the straight line x = 1 in two 
points (r o ,0o) and (ri,0i) with O G (— tt/2,0), 0i G (0,7r/2), r k cosd k = 1 for 

= 0,1, and 7'o > 1 and 6*i must satisfy the equations ©■ 

In [T3] it is proved that the discontinuous differential equation <JSJ> has three limit 
cycles (i) and that the their points {ro,9o) and (ri,9i) are approximately for the 
inner limit cycle of Figure [T] 



(14) r = 1.013330663139.., 6 = 0.162383740477.., 9 l = 0.5541676264624..; 



3. Application 



(13) 






12 



J. LLIBRE, D.D. NOVAES AND M.A. TEIXEIRA 



for the middle limit cycle of Figure [T] 

(15) r = 1.003945075086.., 6 = -0.088680876377.., 6i = 0.768002346543..; 
for the external limit cycle of Figure Q] 

(16) r = 1.111870463116.., 9 = -0.452434880837.., Q x = 1.034197922817... 

It is easy to check that (fT4")) . (|15p and p^]) satisfies the two equations @. Hence, 
Theorem [2] is proved. 

Appendix A: Basic results on the Brouwer degree 

In this appendix we present the existence and uniqueness result from the degree 
theory in finite dimensional spaces. We follow the Browder's paper [3], where are 
formalized the properties of the classical Brouwer degree. 

Theorem 9. Let X = W l = Y for a given positive integer n. For bounded open 
subsets V of X , consider continuous mappings f : V —> Y , and points yo in Y 
such that j/o does not lie in f(dV) (as usual dV denotes the boundary ofV). Then 
to each such triple {f,V,yo), there corresponds an integer d(f,V,yo) having the 
following three properties. 

(i) If d(f,V,yo) 7^ 0, then y S f(V). If fo is the identity map of X onto Y, 
then for every bounded open set V and yo £ V , we have 

d(f \ v ,V,y )=±l. 

(ii) (Additivity) If f : V Y is a continuous map with V a bounded open set 
in X, and Vi and V2 are a pair of disjoint open subsets of V such that 

yo i f(V\(Vi u v*)), 

then, 

d (fo, V, y ) = d (f , V!,yo) + d (/„, Vi,yo) ■ 

(iii) (Invariance under homotopy) Let V be a bounded open set in X , and con- 
sider a continuous homotopy {ft : < t < 1} of maps of V in to Y . Let 
{yt : < t < 1 } be a continuous curve in Y such that yt ft(dV) for any 
t£ [0,1]. Then d(ft,V,yt) is constant in t on [0,1]. 

Theorem 10. The degree function d(f,V,yo) is uniquely determined by the three 
conditions of Theorem^ 

For the proofs of Theorems [9] and [10] see [3] . 

Appendix B: Basic results on averaging theory 

In this appendix we present the basic result from the averaging theory that we 
shall need for proving the main results of this paper. For a general introduction to 
averaging theory see for instance the book of Sanders and Verhulst [19] . 



Theorem 11. We consider the following differential system 
(17) x'(t) = eF(t,x) + e 2 R(t,x,e), 
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where F : K x D — > W 1 and R : M x U x (—£/,£/) — » E™ are continuous functions, 
T -periodic in the first variable and D is an open subset of W 1 . We define the 
averaged function f : D — > W 1 as 

(18) f(x)= t F(s,x)ds, 

Jo 

and assume that 

(i) F and R are locally Lipschitz with respect to x; 

(ii) for a £ D with /(a) = 0, there exist a neighborhood V of a such that 
f{z) ^ for all z e V\{a} and d B (f,V,0) ^ 0. 

Then, for \e\ > sufficiently small, there exist a T -periodic solution x(t,e) of the 
system (1171) such that x(0,e) — > a as e — > 0. 

Theorem [TTJ for studying the periodic orbits of continuous differential systems 
has weaker hypotheses than the classical result for studying the periodic orbits of 
smooth differential systems, see for instance Theorem 11.5 of Verhulst [23], where 
instead of (i) is assumed that 

(j) F, R, D X F, D X F and D X R are defined, continuous and bounded by a con- 
stant M (independent of e) in [0, oo) x D, —£f < e < £/; 

and instead of (ii) it is required that 

(jj) for a 6 D with f{a) we have that J/ (a) ^ 0, where J/ (a) is the Jacobian 
matrix of the function / at the point a. 

For a proof of Theorem [TT] see [JJ section 3. 
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